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Abstract 

We discuss BSDE with drivers containing nonlinearities of the type 
p(y)|z| and p(y)|z| 2 with p a polynomial of any degree. Sufficient con- 
ditions are given for existence and uniqueness of solutions as well as 
comparison results. We then connect the results to the Markovian 
FBSDE setting, discussing applications in the theory of PDE pertur- 
bation and stating a result concerning a large deviations principle for 
the first component of the solution to the BSDE. 

2010 AMS subject classifications: Primary: 60H10; Secondary: 60J60, 
60F10, 35B20. 
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1 Introduction 



In the past 10 years there has been an explosion of publications in the field of 
backward stochastic differential equations (BSDE) mainly due to their signif- 
icance in optimization problems in stochastic control theory. BSDE provide 
a full stochastic approach to control problems which are usually found in 
the Hamilton- Jacobi-Bellman formalism. Moreover, in a Markovian setting 
they provide a connection to certain classes of parabolic PDE via the gener- 
alized Feynman-Kac formula. An application of this falls into the realm of 
numerical analysis, where it may be better to employ a probabilistic solver 
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for a PDE than a standard deterministic one (because of dimensionality for 
instance). BSDE with driver that grow quadratically in the control vari- 
able (introduced in [KobOOj and referred to as gqBSDE) are of particular 
importance as they appear naturally in many problems from mathematical 
finance, for example in the context of utility optimization with exponential 
or power utility functions in incomplete markets. 

In this work we study an extended class of qgBSDE where we allow 
for cross nonlinearities in the driver function. The canonical setting for 
qgBSDE, apart from a bounded terminal condition, is a driver satisfying a 
growth condition of the type 

\f(t,y,z)\ <C(l + |2/| + |z| 2 ) with C>0. 

Here we extend, for the first time, the framework to a growth condition of 
the type 

\f(t,y,z)\ <C(l + \y\ + (l + \y\ k )\z\ + (l + \y\ k )\z\ 2 ), k E N, C > 0, 

which coupled with the corresponding modulus of continuity condition allows 
us to obtain existence, uniqueness and comparison results. In a second step 
we associate with a BSDE of this type a standard forward diffusion process 
(these systems are called forward-backward SDE and denoted as FBSDE) 
and make the link to the PDE framework via a nonlinear Feynman-Kac 
formula. 

The results we present in this work concerning this class of BSDE serve 
as a stepping stone for future research, for instance, under certain condi- 
tions this type of BSDE can be seen as cash sub-additive risk measures as 
described in Section 8 of [EKR09J. Moreover, they also appear in problems 
of optimal investment and consumption when payments of taxes take place 
(see jEFKYllj ). 

In the PDE framework this type of nonlinearities are also of importance. 
These improvements on existence and uniqueness results for BSDE include, 
as particular cases in the PDE framework, results for the Burgers' and gen- 
eralized Burgers' equations. For example, |CS09| ICS10] consider a problem 
in fluid mechanics that links the Navier-Stokes equation for incompressible 
fluid flow and FBSDE with a driver of the type f(y, z) = yz. But they only 
study that specific FBSDE. 

Another example is the equation ut = Au + u\V x u\ 2 for x £ D, t > 
and where D is a smooth bounded domain of W 1 . This PDE relates to an 
FBSDE with a driver of the type f(y,z) = y\z\ 2 that clearly falls in our 
setting. This type of equation belongs to a class of equations used in a wide 
range of applications: in geometry, as a tool to build a harmonic map that is 
homotopy equivalent to another given map (see e.g. |Str96| ). in the theory of 
ferromagnetic materials (the Landau-Lifschitz-Gilbert equation or in models 
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of magnetostriction) and the theory of liquid crystals (see [BDPvH02] and 
references therein). 

The third message of this paper is concerned with the theory of PDE 
perturbation and the same concept for FBSDE. A rough intuition is as fol- 
lows: given a first-order transport PDE for which one is uncertain wether it 
has a solution or one is unable to ascertain its regularity, then one usually 
adds a vanishing regularizing term (typically a vanishing Laplacian) trans- 
forming the first order PDE into a second order parabolic PDE which has 
better properties (or so one hopes). Then as the regularizing term vanishes 
one aims at proving that part of those better properties carry to the limit. 

We discuss the PDE equivalent to a perturbation technique in the FBSDE 
framework. Here we look at a FBSDE system where the forward diffusion 
equation is perturbed by a ^/eadW term and investigate what happens to the 
FBSDE solution as e vanishes. The scope of the results we are able to obtain 
depend on the type of perturbation we work with. With suitable assumptions 
we are able to go indeed far and even show a large deviations principle 
(LDP) for the first component of the BSDE's solution as the SDE's diffusion 
coefficient vanishes (LDPs for SDEs are well known, see e.g. |FW98j . [DZ09J 
or [FK06J). In a subsequent iteration, we give two particular applications 
of the theory we have just developed to emphasize the results and their 
limitations. The relevance of this type of results in finance is related to the 
work of |SS11| . where the authors look at smile asymptotics derived from a 
model using FBSDE as risk measures. Although the authors start with a 
FBSDE formulation, they quickly change to a PDE formulation in order to 
obtain the said asymptotics. The results we propose in this work may prove 
useful to show such asymtpotics without the need of such PDE results. 

Concerning the techniques we use to show the large deviations principle, 
based on the works of |DR07| and [RaiQ6] , they differ quite substantially from 
the techniques commonly used to prove such a result for diffusion processes. 
In the FBSDE case one is able to directly use the contraction principle 
while in the SDE framework, when the diffusion coefficient depends on the 
state process, one is not. The reason is that the SDE's solution process 
when interpreted in the space of continuous functions, although measurable, 
needs not be continuous in that space. A rigorous derivation of an LDP for 
SDEs can be found for example in [DZ09] or |FK 06j. In our case we are 
able to take advantage of the Markovian framework in which the FBSDE is 
immersed and its inherent PDE characterization. This lightens considerably 
the complexity of the proof and, under certain conditions, allows a direct 
use of the contraction principle. 

The work is organized as follows: in Section 2 we introduce the nota- 
tion we will work with and recall for completeness some results on bounded 
mean oscillation (BMO) martingales. Section 3 contains the main theorem 
concerning existence, uniqueness, comparison and estimates for the class of 
BSDE we are introducing. In Section 4 we apply the results to the Markovian 
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FBSDE setting, link the existence of solution of the FBSDE to the viscos- 
ity solution of the corresponding PDE and discuss differentiability of the 
FBSDE with relation to the initial condition of the forward diffusion. We 
close with Section 5 and the results on the perturbation of FBSDE and the 
corresponding PDEs, a large deviations results and two applications. 

2 Preliminaries 

2.1 Spaces and Notation 

Throughout fix T > 0. We work on a canonical Wiener space (£l,J-,¥) 
carrying a d-dimensional Wiener process W = (W 1 ,--- ,W d ) restricted to 
the time interval [0, T], and we denote by T = {3~t)t<stts,T\ its natural filtration 
enlarged in the usual way by the P-zero sets. We shall need the following 
operators, and auxiliary spaces of functions and stochastic processes: let p > 
2, m, n, d G N, Q a probability measure on (fi, F). We use the symbol for 
the expectation with respect to Q, and omit the superscript for the canonical 
measure P. For vectors x = (x 1 , • • • ,x m ) in Euclidean space R m we write 

\x\ = Q^=i( x ') 2 ) 5 > we a l so define 1 = (1, ■ ■ ■ ,1) for the m-dimensional 
vector with all entries equal to 1. We define by sgn(x) : R — > R the sign 
function yielding 1 if x > and — 1 for x < 0. Let Br(xq) denote the ball of 
radius R > centered around the point xq. By 1a we denote the indicator 
function of a set A. 

For a map b : R m — > R d , we denote by Vb its Jacobian matrix whenever 
it exists. To denote the j-th first derivative of the function b(x) with x G R m 
we write X7 Xj b. For a function h{x, y) : R m x M rf — > R we write X7 x h or X7 y h to 
refer to the first derivatives with relation to x and y respectively. A denotes 
the canonical Laplacian operator. We will use the operator <9 t to denote the 
temporal partial derivative. 

We also introduce the following spaces: 

• C^(M m ,R™) the set of fc-times different iable real valued maps defined 
on M m mapping onto R ra with bounded partial derivatives up to order 
k, and C 6 °°(R m ,R n ) = n fe >iC fe fc (R m , R"); We omit the subscript b to 
denote the same set but without the boundedness assumptions. 

• L p (R m ; Q) the space of J^-measurable random variables X : CI — > R m , 

normed by ||X||x,p= L°° the space of bounded random 

variables; 

• 5 p (R m ) the space of all ^-adapted processes {Yt)te[o,T] with values 

in R m normed by \\Y\\ SP = E[sup te[0jT] \Y t \P]p; 5°°(R m ) the space of 
bounded measurable processes; 
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• T-L p (M. m ,Q) the space of all IR m -valued predictable processes (Zt)te[o,T] 
normed by \\Z\\ H p = E<*[( J Q T \Z s \ 2 ds) p/2 ]v ; 

• BMO(Q) or BM02(Q) the space of square integrable martingales $ 
with $o = and satisfying 

\\®\\bmo(Q)= sup || E Q [($) T - (*) T |Jv] \\H 2 < oo, 

T 

where the supremum is taken over all stopping times r G [0, T\. 

• Hbmo the space of M m - valued % p -integrable processes {Zt)te[o,T] f° r 
all p > 2 such that / Q Z s dW s G MO. We define^] ||Z|| Wbmo = 
||J"ZdW|| B MO. 

If there is no ambiguity about the underlying spaces or measures, we also 
omit them as arguments in the function spaces defined above. 

Constants appearing in inequalities of our proofs will for simplicity be 
denoted by C, although they may change from line to line. 

2.2 BMO processes and their properties 

The BMO space is an interesting space of stochastic processes satisfying 

L°° C BMO C n p >iL p . 

For more details on BMO spaces we refer the reader to |Kaz94j . 

In the following lemma we state some properties of BMO martingales we 
will frequently use. 

Lemma 2.1 (Properties of BMO martingales). 

1) Given a BMO martingale M with quadratic variation (M), its stochastic 
exponential £{M)t '■= exp{M^ — \{M)t} has integral 1, and thus the 
measure defined by dQ = £(M)^dP is a probability measure. 

2) For every BMO martingale M , there exists p* > 1 such tha^£(M) G LP . 
Moreover, there exists a constant C p * depending only on p* and the BMO 
norm of M such that for any stopping time r G [0, T] it holds that 

nS(M) P r \F t ] < C p * (£(M) t ) p * (2.1) 

l If Z * W € BMO then one has automatically that Z G H p for all p > 2. This is a 
consequence of the BMO spaces, for more see the BMO results subsection below. 

2 The number p* can be found through the function ^(x) = {l + log 2(x-i) l*^ — ^ 
defined for all 1 < x < oo and verifying linx z ._ ! . 1 + ^(x) — oo and lim a; _ i . 00 ^(x) = 0. In 
other words, if ||Af|| BM o 2 < *(p*), then £(M) £ L v \ see Theorem 3.1 |Kaz94] , 
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3) If \\M\\bmo2 < 1; therefor every stopping time r G [0, T] 

E[ exp{(M) T - <M) r } | JV] < J . (2.2) 

1 II iw II BMO2 

In particular, if Z *W G BMO, then for every p > 1 it holds that 

e[(£ \z s \ 2 d s yj < P l\\z*w\\^ MO => zetf p . 

Moreover, for any p > 1 and an?@ e € (0, 2) 

E expjp^ |Z s | £ ds} <C<oo, (2.3) 
where C depends on p, e and \\Z * W\\^ M q 2 . 

3 Existence, uniqueness and comparison 

We now introduce the class of BSDE we will work with. The main novelty 
is a setting that allows for nonlinearities of the type y k z or y fc |.z| 2 for some 
positive power k > 1. 

3.1 Assumptions 

Take a BSDE satisfying the following dynamics 

Y t = £ + / T /(a, F s , Z s )ds - / T Z s dW s . (3.1) 

We refer to this stochastic equation as BSDE(£,/). We also introduce the 
set of assumption under which we will be working. 

Assumption 3.1. £ is an J 7 ? -measurable uniformly bounded random vari- 
able, i.e. for some M > we have < M; 

Assumption 3.2. / : ft x [0, T] x M x R d -> R is an J 7 -predictable continuous 
function. There exist k G N and positive constants K , 5 such that for all 
(to, t,y,z) fSlx [0, T] x 1 x K d 

l/M,!/,*)! < if (l + |y| + (1 + |y| fc )|z| + (1 + 7 |y| fc )M 2 )- (3.2) 



This result is known as the John-Nirenberg inequality, see Theorem 2.2 in |Kaz94| . 
4 This inequality follows from (12. 2p . since for every e £ (0, 2) and <5 > there exists a 



constant C<s, e such that \z\ E < C<s, e + for all z. 



Assumption 3.3. There exist k G N and a positive constant K such that 
for all t G [0, T], y, y' G R and z, z' € M d # holds P-a.s. that 

\f(;t,y,z)-f(;t,y',z')\ (3.3) 
<K{(1 + (1 + 7 \z\ + Tl^DCIyl*" 1 + \y'\ k - l ){\z\ + \z'\)) \y - y'\ 

+ (1 + \y\ k + \yf + (1 + 7 |y| fc + j\yf)(\z\ + \z'\)) \z-z'\}. 

Remark 3.4. Drivers like f(y,z) = vyz or f(y,z) = y\z\ 2 can be found in 
applications in physics (the first relates to Burger's PDE, see the introduction 
for remarks on the second). Drivers like f(t, y, z) = OfZ + "y\z\ 2 or f(y, z) = 
— (ay + — by~)\z\ 2 with < a < b are found in applications in finance (the 
last one relates to cash subadditive risk measures, see e.g. [EKR09]). 

Remark 3.5. In the above assumptions we write a domination in terms of 
a power \y\ k but since \y\ q < C PiQ + \y\ p for any q < p and some constant C Piq 
it is clear that this includes any polynomial dependence of y up to power k. 

Remark 3.6. In [Tev08], the main theorem holds under his assumption (B) 
that imposes "iV^/j < Const.", a condition which does not allow for yz 
type non-linearities. However, Proposition 1 of [Tev08] partially covers our 
setting by allowing y 2 terms if an extra smallness assumption of the involved 
data is taken, namely that and ||/(-, 0, 0)||l°o are very small. 

For completeness we quote Theorem 2.3 from [KobOOj. This result plays 
a crucial role in proving that a BSDE under Assumption 13. 2l has a solution. 

Proposition 3.7 (Theorem 2.3 of |Kob00| ). Let Assumption \3. 1\ hold and 
assume a continuous J 7 -predictable function f : f2 x [0, T] x M. x M. d — > R 
satisfies for any (u), i)e!lx [0, T] 

\f(u,t,y,z)\ < A + B\y\ + F\z\ 2 , with A,B,F e [0,+oo). (3.4) 

Then the BSDEfeJ) has a solution (Y, Z) e S°° xH 2 . The process Y 

has continuous paths. 

Moreover, there exists a unique minimal solution (Y*,Z*) (respectively a 
unique maximal solution (Y*,Z*)) in the sense that the solution (U,V) of 
BSDE(r], h) where f < h and £ < n (respectively f > h and £ > r/) satisfies 
Y* < U (respectively Y* > U ). 

As is expected from the theory of quadratic BSDE (see |Mor09] . |IDR! 0] 
or |DR11] ). we can conclude that@ Z * W G BMO. 

Lemma 3.8. Under the conditions of Proposition W7A we have Z G TLbmo- 

5 Here we cannot refer to Lemma 3.1 of Mor09] with her Assumption HI instead of 
just H3.4|l . Her extra restriction 7 > b restricts the methods we use in this work. Please 
compare her assumption HI with our Assumption 13.21 



7 



Proof. We only sketch the proof since this kind of argument is known. Take 
t G [0, T]. Assume that HY^oo and H^Hz, 00 are bounded by the same constant 
M and take a constant a satisfying a > 2F. Let r £ [0, T] be a stopping 
time. Applying Ito's formula to the process Yj = expjaYj} between [r, T], 
(|3.4p and the conditional expectation on J- T , we obtain: 



e a5 + 



r 



Y s [aA + aB\Y s \ + (aF - ^-)\Z s \ 2 ]ds 



-aM 



E 



T 



•Ft 



|ZJ 2 ds 



< e 



l + (ai + a5M)r). 



Since a > 2F we have that 2a (a — 2F) > and hence we easily get from 
the definition of the BMO-norm that Z *W £ BMO. Moreover, from the 
calculations we just did, the BMO norm of Z * W is bounded from above by 
a universal constant depending only on ||£||£oo, ||Y||5oo and the constants A, 
B and F. 

We remark as well that as F decreases the upper bound for the %bmo- 
norm of Z also decreases and vice-versa. □ 



3.2 The main results - Abstract BSDE setting 

We now state the main results of this section. We start with an existence 
and uniqueness result. 

Theorem 3.9 (Existence). Let Assumptions COl and \3.S\ hold. Then the 
BSDE (|3.ip has a solution (Y,Z) in S°° x Hbmo- There exists an upper 
bound for Hy^oo independent of^f. 

Moreover, there exists a unique maximal solution (Y, Z) in the sense that 
Y <Y for any other possible solution (Y, Z). 

Theorem 3.10 (Uniqueness). Let Assumptions \3.R \3.H\ and \3.3\ hold and 
take (Y,Z) to be the solution of BSDE (|3.ip . Let HYH^oo < R, \\Z\\-}i BMO < 
for R, R^ > 0, where R is independent of 7 and R^ is not. Define the 
numbed p* = ^ (2K(1 + 2R k )VT + AK(1 + 2jR k )R^) with q* its Holder 

conjugate. Assume further that 7 satisfies (2q* K8jR k ~ 1 ) 2 R? < 1. 
Then BSDE (|3.ip has a unique solution in S°° x Hbmo- 

Remark 3.11. 1) A small clarification is necessary concerning the existence 
of 7 in Theorem \3.10\. since q* and R^ both depend on 7 and there is the 
possibility of a circle argument. One argues as follows, as 7 decreases so 
does the upper bound for \\Z\\u BMO (this is intuitive but see e.g. Lemma 
\3. 8\) . From the definition of ^ in Lemma \2.R the smaller the upper bound 
for \\Z\\n BMO is, the greater p* is and hence the smaller q* is. So, as 7 

6 The function ^ is defined in part 2) of Lemma [2.11 
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decreases so do q* and R? , making it even easier for condition (|3.13p to be 
satisfied. 

2) The proof of Theorem \3.1(h hints that it is possible to obtain the same 
result but under an assumption of the type f(Y,Z) — f(U,Z) < —a\Z\ 2 (l + 
4>(Y, U))(Y — U) with (ft a continuous positive function. This condition would 
then replace the smallness assumption 0/7. We do not explore this direction. 

Finally, we state a comparison result: 

Theorem 3.12 (Comparison). Let i G {1,2}. Define (Y^Z 1 ) G S°° x 
%BMO as the solution of BSDE (|3.ip with terminal condition £ l and driver 
/j. Assume that and /j satisfy Assumptions \3. A \3.2\ and \3.3\ Assume 7, 
satisfies the condition of Theorem \3.1(ft and tha\j\ 

i 1 < e, and h(t, Y t 2 , Zt) < f 2 (t, Y t 2 , Z 2 ) dt dF-a.s., 

Then we have for all t G [0, T] that Y t x < Y t 2 P-o.s. 

Moreover, if either < £ 2 or fi(t,Y t 2 ,Z?) < f 2 (t,Y t 2 ,Z?) in a set of 
positive dt (g) dP -measure then Yq < Y Q 2 . 

Corollary 3.13 (A priori estimate). Under the assumptions of Theorem 
\MBwith~f = 0, it holds for 5f. := /i(- ,Y 2 , Z 2 ) - f 2 {- ,Y 2 , Z 2 ) that for p> 2 

\\y i - y 2 \\% < ci,c p c pr ,^[(\e - e\ + £ \6fs\d S y rq ] 

\\z l -z 2 \\ p HP <c{||*y||^ + ||*y||^ p + ||*y||| l ,E[( J \5f s \d s yy}, (3.5) 

where the constants C vr ^ q i, C p , Cy and C depend only on the constants 
appearing in Assumptions \3. A \3.2\ and \3.3\ and the upper bounds of the BMO 
norms of Z 1 *W and Z 2 * W. Furthermore, the numbers r, q > 1 are related 
as well to the BMO norms of Z 1 *W and Z 2 *W via point 3 of Lemma \2. A 

The following subsections contain the proofs of the above theorems. 



3.2.1 Proof of the Theorem 13.91 - Existence 

As we mentioned earlier, we now need to state and prove a sequence of results 
on a certain family of BSDEs that approximates the BSDE (|3.1j) . 

7 The comparison theorem still holds if instead of fi(t,Y^ , Z%) < /2(i, Vj 2 , Z%) one 
assumes fi(t, Yt , Zl) < f2(t, Yt, Z}). The proof is a straightforward modification of the 
proof we are giving. 
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A "truncation of the identity" function family 

We start by defining a family of smooth functions that truncate the identity 
function, namely 

Definition 3.14. For each n E N let h n : R — )■ R be a continuously differ- 
entiable function with the following properties: 

• |/i n (x)| < n, \h n (x)\ < \x\ and \V x h n (x)\ < 1 for all x E R; 

• h n {x) = x for any x G [— (n — 1), (n — 1)] and h n (x) = n for any x 
outside (—n, n); 

The sequence (fo n )neN converges locally uniformly to the identity function; 
the sequence (V x /i n ) n gN converges to 1 locally uniformly. 

We call the family of functions (/i n )neN a (differentiable) truncation of 
the identitjU. 

A family of truncated drivers and their properties 

With the above defined family of functions we take a driver / satisfying 
Assumptions l3~2l and [3731 and define the sequence (/n)neN through 



/„(*, y, z) = fit, h n (y), z) , for (w, t,y,z)e!lx [0, T] x R x R d , n G N. 



Using the properties of h n (along with \z\ < 1 + \z\ 2 ) and in view of (|3.2p . 
it is clear that there exists a sequence (C n ) rag N such that for any (y, z) £ 



where the latter inequality corresponds to fj3. 2 [) for /„. 

It is also easy to see that for any (y, z), (y', z') G R x R rf and n E N the 
difference \f n (-,y,z) — f n (-,y' , z')\ satisfies condition (|3.3p with constants 
independent of n: just use that |/i n (x)| < \x\ along with the mean value 
theorem on h n and the fact that |V/i n | < 1. 

A family of truncated BSDE and results on them 

With the family of drivers (/ n ) n eN we define a family of approximating BSDE 
obtained by replacing / in (|3.ip by f n so that 



R x R , n E N 



\f n (-,y,z)\ < K(l + C n (l + \z\ 2 )) with limC n = oo 

n— >oo 

\f n (;y,z)\ < K(l + |y| + (1 + \y\ k )\z\ + (1 + 7 \y\ k )\z\ 2 ] 



) 



(3.6) 



(3.7) 




(3.8) 



8 Such a family of functions clearly exists, for an explicit example we point to Subsection 
3.2.1 of [DRllj . 
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The next result states the existence and uniqueness of the solution to (13. 8p . 
The argumentation here follows two steps: in the first we show existence of a 
solution to (|3.8p where we obtain upper bounds for the solution that depend 
on n G N. In a second step, we twist the arguments and obtain upper bounds 
for the norms of 5°° and T~L bmo independent of the truncation height n£l 

Lemma 3.15. Let Assumptions \3.1\ and \ 3. HI hold. For each n G N, BSDE 
(|3.8p has a (unique minimal/ maximal) solution (Y n ,Z n ) in S°° x Hbmo- 

Proof. This follows immediately from Proposition 13 . 7l and Lemma T3 . 81 applied 
to (|3.8p . using that f n satisfies (|3.6p . □ 

The next lemma shows that it is possible to estimate the norms of 
(Y n ,Z n ) solution to (|3.8p in 5°° x Hbmo by universal constants indepen- 
dent of n G N imposing only Assumptions 13.11 and 13.21 We use (|3.7p instead 
of 

Lemma 3.16. Let Assumptions \3.1\ and hold. Let n G N and take 
(Y n ,Z n ) to be the solution of BSDE (|3.8p as given in the previous lemma, 
then 

sup { ||y n || 5 ~ + WWhbmo } < °°- 

In particular sup n ||y n ||s°o < e KT (M + KT), which is independent ofj. 

Proof. Let t G [0, T] and n G N. 

Step 1 - preparation: We start by going back to BSDE (|3~%]t . Let (Y n , Z n ) 
solve BSDE (|3.8p . then we can decompose the growth of driver f n (t, Y™, Z™) 
by using the growth assumption (|3.7p in the following way: 

\f n (t,Y?,Z?)\<K + ftY? + %.Z?, 

where j3 n and b n are two processes valued in R and M d respectively defined 
by: (the number k follows from Assumption 13. 2p 

ft := Ksgn(Y t n ), and b n t := K 1 + l^"l fc + + 7l^ w | fc )l^l Z ?l {{z ^ 0} . 

\ Zj t I 

It is clear that sup neN ,te[o,T\ \\(3 n \\s°° < oo and ((3 n Y n ) = K\Y n \. Moreover, 
b n. Z n = (i + |yn|* + (i + ^|yn|fc)|^»|)|^»| and since Y n G 5°°, Z n G Ubmo 
and (|3.2p we have 

\b n \ <K{l + \Y n \ k + (l + 7 |y n | fc )|Z n |} => b n *W£BMO. 

Hence the probability measure Q n with Radon-Nikodym density dQ n /dP = 
£{b n * W) T is well defined and with relation to which = W. — f Q b™ds 

is a Brownian motion. 
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Step 2 - uniform 5°° bound for Y n : We set e" := exp{ J Q * /3™ds}, which 
satisfies e~ KT < e n < e KT , and use a change of measure to show that an 
upper bound for ||y n ||s°° can be obtained independently of n (and 7). 

We take (Y n ,Z n ) as the solution of BSDE and with the help of 

processes (3 n , b n and e n defined above we write the BSDE for e™Y™ under 
<Q n via Ito's formula. We have then 

e nyn<^_ [ T e n Z n dWs+ F e »[-ftY? + \ f n (s,Y?, Z?)\]d8 
Jt Jt 

< e n T i - F e?Z?dW a + F e?HW + K + ^Y s n + b n s ■ Z n s ]ds 

Jt Jt 

< e KT M - F e n s Z™dWT + e KT K(T - t). (3.9) 



where we used the properties of the process e n , the measure change to Q n 
and that £ G L°°. We obtain by taking Q"-conditional expectation that for 
some positive constant C and for all t £ [0, T] 

Y t n <e KT (M + KT) =>■ sup \\Y n \\ S oo < e KT (M + KT). 

neN 

where the last conclusion follows from using similar arguments to obtain a 
lower bound for Y n . 

Step 3 - uniform %bmo bound for Z n : The argument uses that we already 
know that sup nGN ||y n ||5°o < 00. This yields the existence of n* E N such 
that for any n G N with n > n* 

f n (t, Y n , Z n ) = f n *(t, Y n , Z n ) = f(t, Y n , Z n ). (3.10) 

The arguments used to prove Lemma 13.81 imply that the upper bound for 
the BMO norm of Z * W depends only on the problem data and an upper 
bound for ||y||sao. Hence we get here that sup ngN ||Z n ||% SMO < 00. □ 



The proof of the existence theorem 

We are now ready to prove the existence theorem. Lemma 13.161 yields a se- 
quence (Y n , Z n ) n< z.^] of processes belonging to 5°° x Hbmo that solve BSDE 
(|3.8p for each n G N. Moreover, the processes' respective norms are bounded 
uniformly in n. We need only to conclude that a limit for the sequence 
(y n ,Z n ) ng N exists and that the said limit solves BSDE (|3.ip under As- 
sumption 13.21 

Proof of Theorem VJ.ih In view of (|3.10p it is clear that for each n > n* pairs 
(Y n ,Z n ) solve both the BSDE <J3T8]) and (pTT]l . And hence we obtain the 
"existence of solution" result. 

Moreover, since for each n > n* there exists a maximal solution to (|3.8p . 
we have the existence of a (infinite amount of) maximal solutions to (|3.1j) . 
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It is clear that the maximal solution is unique, in the sense given in the 
theorem's statement. If (Y, Z) and (Y,Z) are two maximal solutions then 
we have that Y < Y and also Y < Y. This implies also the uniqueness of a 
maximal solution of BSDE (13.11). □ 



3.2.2 Proof of Theorem 13.101 - Uniqueness 

We now present the proof of the uniqueness result. As usual to prove ex- 
istence, one needs only bounds on the growth of the involved functions but 
to prove uniqueness one needs to control the modulus of continuity of the 
involved functions, hence for uniqueness one additionally needs Assump- 
tion 13.31 

Unfortunately, the arguments we present here do not allow for general 
choices of 7 (from Assumption 13.21 or I3.3|) . Here 7 has to be small enough. 
This smallness enters in play in f|3. 14[) below. In the later Subsection 13.31 we 
discuss a particular situation where it is possible to solve the BSDE for any 
value of 7. 

Proof of Theorem \3.1(K Let t € [0, T] and assume that (Y,Z) and (U, V), 
belonging to the space 5°° x HbmOi are two solutions to the BSDE (|3.ip . 
The 5°°-norms of Y, U are universally bounded by a constant, say R, inde- 
pendently of 7 (see Lemma 13. 16[) : the Hbmo- 110 ™ of Z, U are bounded by 
a constant R^ that depends on 7 and decreases as 7 decreases (see the proof 
of Lemma 13. 8|) . 

Now, define 5Y = Y — U, 5Z = Z — V and the processes 



_ f(t,Y t ,Z t )-f(t,Ut,Z t ) 
1 ' Yt-U t 



l{y t -ET t #o}. et := exp j J r s ds|, (3.11) 



~ nt,u t ,z t )-f(t,ut,v t ) /r7 

and h = | Zf _ prp y Zt ~ Vtl^Zt-Vt&o}- (3-12) 

Concerning the above defined processes, we remark that for b we have from 
(|3.3p (and the triangular inequality of the ^bmO" 11 ™) that 

|S| < K{\ + 2R k + (1 + 2 7 R k ){\Z\ + \V\)}, 
\\b\\H BMO <C := 2K{(1 + 2R k )VT+ 2(1 + 2 7 R k )R^}. 

Since Z,V £ Hbmo, it is easy to see that the process b*W € BMO and 
hence the probability measure Q with Radon-Nykodim density dQ/dP = 

£(b * W) is well defined and = W. — f Q b s ds is a Q-Brownian motion. 

Moreover, from Lemma [2. II and defining p* = ^ _1 (C) with q* its Holder 
conjugate we have that £ (b * W) E LP . 

In view of fl3J3} we have that |r| < K(l + AR k ~ 1 \Z\ +87i? fc " 1 |Z| 2 ). Since 
Z £ Hbmo-, the properties of BMO martingales, namely (|2.3p . combined 
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with the behavior of R? (described in the proof of Lemma 13. 8|) imply that 
it is possible to find a 7 small enough such that|f|: 

(V**^" 1 )* < -L < ^ . (3.13) 

R~ \\Z * W \\bmo 

In view of (|2.2p and the strictness of (|3.13p it is always possible to find a 
small enough e > such that 

if G L 1+£ . (3.14) 

It is also clear that e has continuous paths. 

Step 1 - Uniqueness of the solution's first component: With the help of 
the processes b and e defined above we proceed as in step 2 of the proof of 
Theorem 13.91 (we skip some details and point the reader to (|3.9p ). Using 
Ito's formula, we write a BSDE for (ifSYt) under the Q-measure: 

e t 5Y t + ^ e s 5Z s dwf 

= ° + [ T U-TsSY s + f(s, Y s , Z s ) - f(s, U s , Z s )]ds = (3.15) 

We obtain it SYt = Q- and P-a.s. for all t G [0, T] once we justify that 
the stochastic integral f e s SZ s dWg is indeed a true Q-mar ting ale. For this 
we have to essentially prove the square integrability of the process, i.e. the 
finiteness of 



E Q [ [ \e s \ 2 \5Z s \ 2 ds] <E\£(b*W)\e T \ 2 [ \6Z s \ 2 ds 
Jo L j . 



< E 



£V>* FT) )' p " E\\e T \ 2q * ( f T \5Z s \ 2 ds^ q 



Combining now with the fact that £{b * W) G LP , that bZ s G Hbmo an d 
hence 5Z G T~L P for all p > 2 and finally (|3.14p . we can obtain the sought 
conclusion. 

At this point it is clear that etSYt = for every t G [0, T] P-a.s. We now 
want to conclude that we also have 5Y = 0. The continuity of e (and 8Y) 
yield that there exists a set A £ Q satisfying ¥[A] =0 such that 

for all t G [0, T], 00 G Vt\A we have e t (u)) > 0. 

Now, given the positivity of (et)te[o,T] we can conclude that Yt = Ut for any 
t G [0, T] P-a.s. 



9 In short |r| < -y\Z\ 2 and we are faced with the question of the integrability of e. From 
the John-Nirenberg inequality in (|2.2|l we are only able to ensure the integrability of e 
when the constant 7 is sufficiently small. In general the quadratic variation of a BMO 
process Z * W is not exponentially integrable. 
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Step 2 - Uniqueness of the solution's second component: We are missing 
only the uniqueness proof for the control component of the BSDE. Here we 
return to (|3.15p and take advantage of the fact that we already know that 
SY = 0. We have, by using ito's Isometry, that E®[Jq \e s \ 2 \5Z s \ 2 ds] = 0, 
from here using the positivity of (et)t^\o ; T] ( see the above step), we conclude 
that \5Zt\ = \Zt — Vt\ =0 dt ® P-a.s. This concludes the proof. □ 



3.2.3 Proof of Theorem 13.121 and Corollary 13.131 - Comparison 

We now prove Theorem 13.121 and Corollary 13.131 The first parts of these 
proofs are very similar to that of Theorem 13.101 and so we do not give the 
full details. 

Proof of Theorem\3AE Let t 6 [0, T]. Define 5Y = Y 1 — Y 2 , 5Z = Z 1 - Z 2 , 
5£ = e- e and 8f. = fi(-,Y 2 , Z 2 ) - / 2 (-, Y 2 , Z 2 ). 

We define as well, processes similar to those of (|3.1ip and (|3.12p . namely: 

Tt ■= Y 1 -Y 2 h*? -*?¥*>}> et := exp \ J r * ds /' 

(3.16) 

, r _ h(t, Y t-, z t) ~ Y t-, z t) 1 7 \ 7 i\a ("?17^ 

and °t - | Z i _ Z 2|2 ( Z * ~ Z t ) 1 {[Z t 1 -Z||^0}- l 3 - 17 ) 

And as before, with b * W we define a new probability measure Q with 
Radon-Nikodym density dQ/dP = £(b * W) and W® = W. - J Q b s ds is a 
Q-Brownian motion. 

From the BSDE for SY, and Ito's formula we can write a BSDE for the 
process etb~Y t as 

e t 5Y t = e T 5£ - [ e s 5Z s dWf + [ e s 5f s ds 
Jt Jt 

= ¥^\e T 5i + J e s Sf s ds\^. (3.18) 

Using that et > and the theorem's assumptions, <5£ < and 5ft < 0, 
we conclude that etSYt < and hence that for all t £ [0, T) it holds that 
5Y t = - Y t 2 < Q-a.s. and hence also P-a.s. 

We see also that at t = if 6£ < or if 5f < in a set of positive 
dt (g) dP-measure, then the inequality is strict and we conclude that 5Y$ = 
Yq — Yq < 0, which proves the result. □ 

We finish this subsection with the proof of Corollary 13.131 This proof is 
close to that of Lemma 3.2 in [IDR10J but with a different argumentation in 
what the linearization trick is concerned. 
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Proof of Corollarv \3.13[ 

Step 1 - The estimate for Y 1 — Y 2 : To prove this corollary we stop at the 
point of inequality (|3.18p in the proof of Theorem 13.121 and continue it from 
there but in a different fashion. Recall that e and T were defined in (|3.16p 
and (|3.17p . and from the corollary's statement, it is assumed that 7 = 0, this 
means that (|3.14p now holds without any restriction on p. 

We remark as well that for any < t < s < T 



exp 



f r dr} < A T := exp { J +4R* _1 |Z a 1 |)d*}. 



Since Z 1 G HbmOi t ne random variable At is integrable in view of (12. 3p . 
To shorten the notation, we define St := S{b * W)t as the density of the 
probability measure Q. We have then 



8Y t < E K 



< E % 



e T (e t )- l m + (e t )- 1 / e s \Sf 8 \ds 



T 



A T [m+ \5f 8 \ds 



ft 



From the above inequality, it follows that 



5Y t < S t - 1 E ¥ \S T A T (|^| + / \6f,\ds 



<S t - 1 E ¥ [S T ^E F \AT(\5Z\+ / \6f s \ds 



ft 



< C r ,W 



A r T m + / \Sfs\ds 



T 



where the last lines follow from b * W G BMO combined with (|2.ip and 
Holder's inequality (with 1/r + 1/r' = 1) for r' > 1 with S T G L r '(P). 

Let now p > 1 and apply Doob's inequality (and concave Jensen's in- 
equality) to conclude that 



E[ sup \5Y t \P) < C p r ,C p E A p T r (m + 
te[o,T\ 1 v 



\Sf s \ds 



pr 



< C P ,CpCp r ' t q'E 



m + / \5f s \ds) 



\ prq 



where we applied once again Holder's inequality with exponents 1/q + l/q' = 
1 and version of for q' > 1 with A% G L«'(P). 

Step 2 - The estimate for Z 1 — Z 2 : We follow the notation introduced 
above. The estimate here is obtained using standard BSDE techniques, so 
we shorten the proof a little. We first remark that 



\M;Y l ,Z l )-M;Y. 2 ,Z 2 )\ <U\5Y.\+L Z \5Z.\ 
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where, using fl3J3]) with 7 = and that {{Y^f < R for i 6 {1, 2}, 

L y := K(l + 2R k ~ l (\Z l \ + \Z 2 \) and L* := K{\ + 2R k + |Z X | + \Z 2 \). 

Applying Ito's formula applied to |<5y| 2 over the interval [0, T], using the 
"Lipschitz" assumptions of /, Young's inequality (with parameter 2) and 
observe that with 5f defined as 5f. = Y 2 , Z 2 ) - / 2 (-, Y 2 , Z 2 ), we obtain 

|<5y | 2 + [ T 26Y s 5Z s dW s 
J 

= |^| 2 + / [25Y s {f x {s,Y},Zl) - f 2 (s,Y 2 ,Z 2 )) -\5Z s \ 2 }ds 
Jo 

< \5C\ 2 + [ [2(L» + (^) 2 )|<5y s | 2 + 2|<jy fl | \6f s \ - l\6Z s \ 2 ]ds 
Jo 1 

1 c r / \ 2pi - 

=► < C{py|| 2p 2p + ||5Y||g p E[( ^ [L v s + (L z s ) 2 ]ds) ] 2 

\8f s \d s y P ] h +e\\5Y\\% p + -\\8Z\\ 2 ^ p }, 



+ \\SYf s2p E 



where the last line follows from: reordering the terms, taking absolute values, 
exponentiating to power p and taking expectation, then use Ito's Isometry, 
the fact that 5Y £ S°° and Young's inequality for some e > 0. Since 
Z l 6 %bmo f° r * £ {1)2} the processes L V ,L Z G Hbmo, now with the 
appropriate choice of e, estimate (|3.5p follows. □ 

3.3 A remark on Theorems 13.91 and 13.101 

Although we managed to show that a solution to BSDE (|3.ip exists under 
Assumption 13.21 for any 7 > 0, we were only able to show a comparison 
result when 7 was small enough (see (|3.14p ). This is due to a limitation of 
the mathematical tools available, namely that a small 7 ensured that the T's 
in (|3.1ip and (|3.16p had the required integr ability properties. 

But there are situations where one is able to show uniqueness for any 
value of 7. We next discuss about the possibility of linearizing a driver, 
e.g. f(y,z) = g(y)\z\ 2 with a convenient g function, through an invertible 
transformation. The idea is similar to the usual quadratic driver case, f{z) = 
\z\ 2 /2 and the transformation Pt = exp y. 



A linearization transformation 

Linearization for standard quadratic BSDE is well understood: through the 
transformation Pt = exp{7y} and Qt = ZfPt, the BSDE 

Y t = Z + f T l\ Z s\ 2ds - [ T ZsdWs 
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becomes 

Pt = e^- [ Q s dW s => Y t = -logE[e*|.F t ]. 
Jt 7 

It is possible to use the same type of arguments and find linearizing transfor- 
mations for BSDE with drivers of the form f(y, z) = g(y)\z\ 2 for an integrable 
function g, by solving the ODE0 

g(vW(y) - \&'{v) = o, $(o) = o, <&'(o) = 1. 

The above ODE is solved by 

<f>'(y) = exp | J 2g(x)dx^ and <J?(y) = ^ $'(x)dx. 

Since 5 is integrable by assumption, we have <&'(y) 7^ for any y 6 R so 
that $ is invertible and hence the unique solution of the transformed BSDE 
yields the unique solution of the original BSDE. 

We also point out that this type of transformation works well in the one 
dimensional setting. In the multidimensional case one quickly finds problems, 
see the Remark 14.91 below. 



4 FBSDEs and PDEs 

In this section we discuss the framework of Markovian FBSDE and their 
connection to PDEs. The results here open way to the later results on PDE 
perturbation and large deviations. 

Here we particularize the previous results to the framework of decoupled 
Forward-Backward SDE (FBSDE), where we assume that the randomness of 
the terminal condition and driver have origin in a diffusion process. Namely, 
we define an FBSDE, with solution (X,Y,Z), as the following system of 
stochastic equations: let m £ N and take (t, x) £ [0, T] xM™ and for s £ [t, T] 

Xl' x = x + j b(r,X t r ' x )dr + j a^X^dWr, (4.1) 
Y^=g{X^)+ f f(r,X^,Y r ^,Z^)dr- f Zj> x dW r . (4.2) 

J s J s 

The involved functions b, a, g and / satisfy: 

Assumption 4.1. b : [0,T] x R m -> R m and a : [0,T] x R m -> M. mxd 
are continuous functions for which there exists a constant K > such that 
su Pte[0,T]{ IK^)0)l + I , 0) I } < K. Moreover b and a satisfy standard Lip- 
schitz conditions in the spatial variables. 



10 Apply Ito's formula to &(Y t ) and organize the terms inside the Lebesgue integral in 
a convenient way. Note that we already know that Y £ S°° and can easily verify that 



f <&' (Y)ZdW is a true martingale. 
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Assumption 4.2. g : R m — >■ R and f : [0, T] x R m xftxR" 1 -*! are conim- 
uous functions, g satisfies a standard Lipschitz condition and is uniformly 
bounded by a constant M . Let k G N and a positive constant K such that 
for all (t, x, y, z) G [0, T] x R m x R x R d 

\f(t,x,y,z)\ <K(l + \y\ + (l + \y\ k )\z\ + \z\ 2 ). (4.3) 

Let t G [0, T], x, x' G R m , y, y' G R and z, z' G R d it holds that 

\f(t,x,y,z)-f(t,x',y',z')\ 

< K{(l + \y\ + \y'\ + \y\ k \z\ + \y'\ k \z'\ + \z\ 2 + \z'\ 2 )\x 
+ (l + (|y| fc - 1 + |yf- 1 )(|z| + |z'|))|y-y'| 
+ {l + \y\ k + \y'\ k + (\z\ + \z'\))\z-z'\} 

4.1 Particularizing the main results - Decoupled FBSDE 

We now state a particularization of the main theorems in the previous section 
to the framework of FBSDE. 

Proposition 4.3. Take (t, x) G [0, T] x R m and let Assumptions ^. 1\ and \4-S\ 
hold. Then there exists a unique triple {X t)X ,Y t,x , Z t,x ) G S p x 5°° x Hbmo 
for any p>2 that solves FBSDE f[47T]) - fl472]) . Moreover, 

sup H^'ls- + \\Z^ x \\n BM o < oo. (4.5) 

(t,i)6[0,T]xI m 

Furthermore, Y t,x and Z t,x are Markovian: there exist two Borel-measurable 
functions u : [0, T) x R m — x R and v : [0, T] x R m — >■ R d such that 

Y s f ' x = u{s,X t s ' x ), s G [t,T], P-a.8 Z l f = v(s, X l f), &s®R-a.s. 

Proof. Existence and uniqueness follow from standard SDE results and the 
results of the previous section. 

Let (t,x) G [0,T] x R m . The uniform bounds of the and Hbmo- 
norms, with relation to t and x, follow from the uniform bounds of the growth 
assumption in the driver / and the arguments used to prove Theorem 13.91 

The constant K appearing in the domination (|4.3p is independent of x 
and t and it is shown in the proof of Theorem 13 . 9 1 that the uniform bounds for 
the norms of Y and Z depend only on the constants K and M appearing in 
Assumptions 13. ll and 13.21 which are essentially the same constants appearing 
in Assumption 14.31 

The Markov property of (Y, Z) is also trivial. It follows from arguments 
similar to those in Lemma 4.1 in |EKPQ97,] or Theorem 4.1.1 in |DR11| . □ 

The next result states the continuous dependence of the solution in the 
Euclidean parameter x. 



(4.4) 

-x'\ 
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Corollary 4.4. Take t £ [0, T] and x,x' £ M. m and let Assumptions \4 . 1] and 
4-^ hold. Then for any p > 2 there exists a constant C > ( independent of 
x and x' ) such that 



E[ sup \Y*' X - Y*' x '\ p ] + E ( / \Z l f - Zl' x '\ 2 ds ) 2 <C\x- x'\ p 



\Y*' x -Yy\ p ]+E ( [ T 

se[t,T] 

For < t < s < T the mapping M m 3 x i— )■ l^*' 1 admits a continuous 
modification where almost all sample paths are a-Holder continuous in W 71 
for any a G (0,1). For G [0,T] x M m i/ie mapping s \— > Ys' x (lo) is 

continuous for P-a.s. a; G 17. 

Proof. For simplicity we take t = 0. The bounds (|4.5p play a crucial role 
because we can now apply the results of Lemma 13.131 where the involved 
constants C are independent of the parameters x and x' . Then, for p > 2 
we have (recall that H^Hsoo is uniformly bounded in x) 

\\v x — v x ' \\p 

w 1 1 \\sp 



< CE 



(\g(X%) - g{X$)\ 
+ £ \f(s, X X ,Y X , Z x ) - f(s, X X ',Y X , Z x )\ds) 

< C I \ \X X — X X ||^pqr + 



+ E 

' - 'te[o,T] 



sup \Xf-Xf\ I [l + |Zf| + |^| 2 ]d S 



Using \z\ < 1 + \z\ 2 along with the fact that ||-^ x ||wsmqJ, s uniformly bounded 
in x, the result follows from standard results of SDH 11 !, namely that under 
Assumption 14. II we have for any x,x' £ M m and p > 2 that 

\\X X - X X '\\ P SP = E[ sup \Xf - X X '\P] <C p \x- x'\ p , 
te[o,T] 

with C p a constant independent of x and x'. The continuity of the map in 
t and x follows from the properties of the solution of the BSDE and Kol- 
mogorov's continuity criterion for p large enough (see Theorem 2 in [SchOOj). 
So we can deduce the result concerning \\Z X - Z x from ([33]). □ 

The previous result hints that it is perhaps possible to make sense of the 
variation of Y x , i.e. that x t— >■ Y x is differentiable (in some sense). 

Using a particular form of the results in |BC08] we can conclude the 
existence of all partial derivatives (in x £ M m ) of (Y X ,Z X ). What is then 
left to argue is the existence of the total derivative of Y. We will work under 
the following assumption 



'See for instance Theorem 1.2.5 of [DRllj 
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Assumption 4.5. Let Assumptions ]^, i] and \4.2^ hold. b, a and g have con- 
tinuous, uniformly bounded first order spatial derivatives, f is differentiate 
with continuous derivatives satisfying for some K > for all (t, x, y, z) G 
[0, T] x W n xlx m d 

\V x f(t,x,y,z)\ < K{1 + |y| + \y\ k \z\ + |z| 2 ), 
\V y f(t,x,y,z)\ <K{l + \y\ k - l \z\), 
\V z f(t,x,y,z)\<K(l + \y\ k + \z\). 

Let t E [0,T], x G R m , 9* = (F, Y x , Z x ) be the solution of (l4TT1)-(l42l). 
and define, at least formally, the dynamical 



(4.6) 



V x Xf = I d + f V x b(s, X x )V x X x ds + f V x a(s, X x s )V x X x dW t 
Jo Jo 

V x Y t x = (V x g)(X%)V x X%- J V x Z x dW s (4.7) 

+ [ <(v/)( s , e x s ), (v x x x ,v x y x , v x z x ))ds. 



Proposition 4.6. Let x G W" 1 , p > 2 anc? /e£ Assumption ]!^. 5] hold. Then the 
mappingW 1 -> SP(M m ) xSP(M) x?P(M d ) ; x (-)• (X X ,Y X ,Z X ) is differentiate 
in the norm topology and the derivative is the unique solution of FBSDE 
gSl)-g2| in S p xS p x ?P /or any p>2. 
In particular, for x', x G M m we /tawe 

jim {iiv.y^ - v.y^'n^ + 1| v x z*> x - v x z^ x '\\ p HP } = o. 

Remark 4.7. As m Theorem 2.2 of !AIdB.07f (or Section 3.1 of [DEI If ) it 

is possible to show that there exists a function £1 x [t,T] x W 71 — > R mapping 
(uj,s,x) i — y Ys ,x such that for almost alluj, Ys' x is continuously differentiate 
in x, i.e. understanding x i— > Y s t,:c (a;) as a function in C 1 (W n ,M)! 

For such a result one would need to ask that Vg satisfies a Lipschitz 
condition as well a Lipschitz like condition for the derivatives of f. Since the 
argumentation for such a result follows from a straightforward adaptation of 
the proofs in IAIdR07l or [DR11], we refrain from giving the details. 

The proof of Proposition 14.61 follows a methodology rather standard for 
this type setting and so it is given in the appendix. 

4.2 Connection to PDEs 

In the classical theory of FBSDE one is able to connect the solution of 
FBSDE (ECS) to tlie solution of a certain PDE. Formally, for < t < 



12 For m £ N, I m represents the m x m-dimensional identity matrix. 
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s < T, xeR m the FBSDE (|4TTjl - (j43|l is related to the PDE: 

d t u(t,x) + Cu(t,x) + f(t,x,u(t,x),(V x ua)(t,x)) =0, (4.8) 

u(T,x) = g{x), 

with the second order differential operator C being the infinitesimal semi- 
group generator of the Markov process (Xt' x ) se [ t x] given by 

m ^ m 

£ ■= y^bj(t,x)d Xi + - ^2 {°°' T )i,j{tiX)dx l x j - 

i=l M=l 

The relation between the FBSDE and the PDE is given by the identities: 

Y*' x = u(s,Xl' x ) and Z?f = (V x u a){s, X l f), (4.9) 

where the second relation holds provided the derivatives of u are appropri- 
ately defined. Indeed, under Assumptions 14 . 1 1 and 14. 2 [ we get from Corollary 
PI that 

\u(t,x) - u(t,x')\ p = \Y*' X - Y l /\ p < C\x - x'\ p , (4.10) 

and hence u(t, x) is a Lipschitz function in its spatial variables. With the 
results from Proposition 14.61 we have also that V x u exists and is continuous 
(in the spatial variable). 

Remark 4.8. Particularizing the forward diffusion by assuming 6 = 0, a = 
y/2 and f(t,y,z) = vyz, one obtains the known Burgers' PDE. 

Remark 4.9 (Multidimensional case of f(y, z) = y\z\ 2 ). Related to a partic- 
ular case of the Theorem we make a small remark concerning the PDE 
representation of such BSDE (when paired with a SDE) in higher dimen- 
sions. There exists a phenomenon of finite-time gradient blow-up in PDEs: 
[CDY92] (also presented as Theorem III. 6. 14 in (Str08]) consider mappings 
u from the closed unit disk in M 2 into the unit sphere in M 3 which satisfy 

d t u = Au + \V x u\ 2 u, u(0,x) = u (x), u(t, -)\ dD 2 = u o\ dD 2- ( 4 -ll) 

They show that for some smooth and bounded boundary condition uq, the 
solution of blows up in finite time, i.e., the maximal existence interval 

[0, T) has a finite T. 

Moreover, in [BDPvH02} for the same PDE 

d t u = Au + \V x u\ 2 u in Q := Q x R + 

u = uq on dQ. x R + and u(0, x) = uq(x) for 

with £1 a bounded domain in M. n with smooth boundary, the authors show that 
uq G C°° exists such that an infinite number of solutions to the PDE exist. 
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4.3 Viscosity solution 



As far as the relation between PDE and FBSDE go, the natural counterpart 
of the FBSDE solution is a PDE solution in the viscosity sense via (|4.9p . 
Observe that we have not mentioned in any way an ellipticity assumption 
for <r, which is consistent with viscosity theory. Although we have a unique 
solution for the FBSDE, the FBSDE is only to provide existence for the 
PDE while uniqueness of the PDE has to be shown via a comparison result 
for PDE. It is unclear whether such a result holds in this framework that 
allows for cross term nonlinearities like yz since this type of nonlinearities 
implies that the Hamiltonian arising from (|4,8p is not a proper function. 
Nonetheless, such a comparison result exists in the probabilistic framework 
(see our Theorem 13. 12[) . 

We do not give many details below, but as expected, one is able to prove 
that the Y component of the solution to the FBSDE (|4.1]) - ()4.2|) is indeed a 
viscosity solution to PDE 



Theorem 4.10. Take (t,x) € [0, T] x M. m and let Assumptions and\J^ 
hold. Then the mapping u : [0, T] x W 71 R defined by u(t,x) := Y^' x is a 
uniformly bounded continuous function of its temporal and spatial variables 
and a viscosity solution of (|4.8p . 



The proof of this theorem is the same as that of Theorem 3.2 in [ Par99| 
except for the application of the comparison theorem in the last step of 
the proof (where we use our Theorem 13. 12[) . We postpone the proof to the 
appendix. 



5 PDE Perturbation and large deviations 

In this section we look at the relation between FBSDE and PDE perturba- 
tion. In a first approach, the form of the perturbation is quite amenable to 
a in-depth analysis including the possibility of showing a large deviations 
principle for the solution of the FBSDE as the perturbation vanishes. The 
second step looks at another setting of PDE perturbation and only partial 
results, within the framework of the previous sections, are obtained. We then 
stretch the theory to recover some other known results on a certain perturbed 
transport PDE and the Navier-Stokes equations for incompressible fluid flow. 

In this section we establish a large deviations principle for the FBSDE 
(|4.ip . (|4.2p . This type of result has its relevance in the theory of perturbation 
of PDEs. We discuss two types of perturbations, while the first leads to 
classical results (and relates to the theory of e-vanishing viscosity solutions) 
the second shows the limitations of the techniques used. 
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5.1 Canonical perturbation in FBSDE framework 

Throughout let x G R m , < t < s < T and e G (0,1]. Moreover, we 
denote by C a constant that may change from line to line but it is always 
independent of the parameters e, x or i. 

Throughout we work with Assumptions 14.11 and 14.21 and from the context 
it is clear that the constants involved in the growth bounds are independent 
of e\ At a later stage we will discuss the situation \f(y,z)\ < K\y\ \z\je. 

We work with the family of equations indexed in e and following the 
dynamics 

Xl' x ' £ = x+ b(r,X t r ' x ' £ )dr + ^ o-(r,X$' x ' E )dW r . (5.1) 
Jt Jt 

It is known that as e vanishes, the solution X t,x,£ converges to X t,x, ° (see 
Lemma 15 . 1 1 below) . which solves the deterministic equation 



Xl' x '° = x + ^ S &(r,X*' x ' )dr. 



(5.2) 



Moreover, under some extra assumptions, that the law induced by X t,x,£ 
satisfies a large deviations principle (LDP) (see Definition 15.41 below). 
We want to show that the FBSDE with dynamics 

Y*> x > e = <?(X*f' £ ) - / Z t r > x > E dW r + /(r,A^y^Z^' £ )dr, (5.3) 

J s J s 

converges to the deterministic backward equation 

rT 

Yf> x >° = giX^' ) + / /(r,X^°,y r '- :E ' ,0)dr, (5.4) 

J s 

and that the law induced by Y t,x ' £ satisfies an LDP as well. 
We start with a small result on the properties of X t,x,£ . 

Lemma 5.1. Under Assumption \4- 1\ and for any e G (0, 1] equations (|5.ip 
and (|5.2p have unique solutions, X t,x,£ G S p for any p > 2 and X t,x '° G 
C°([t,T] x M m ,M m ). Moreover, for any p > 2, for some C p > it holds that 

sup ||X*- s ' e |||p < C p {\ + \x\ p ) and ||X*' X ' £ - X^'°||£ p < C p (l + \x\ p )e p . 

£6(0,1] 

Furthermore, lim e _>o X t,x,£ = X t,x, ° a.s. 

Proof. These results are quite standard from either SDE or ODE theory and 
follow from the good properties of the b and a functions. See e.g. |DZ09] , □ 

We next establish properties of (Y t,x ' e , Z t,x ' £ ) and Y t,x, ° . 
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Proposition 5.2. Let Assumptions \4-l\ and \4~J\ hold. Then, a unique so- 
lution (JC*.*» e ,y*» !B » e ,Z t >*' E ) of FBSDE ([511, ([O]) in x 5°° x ft BA/0 for 
any p > 2 exists. Moreover, 

sup sup {||y* ,x,e ||s« + ll^ M,e Hw B Aio} < °°- 

£6(0,1] (t,i)e[0,T]xE™ 

T/ie mapping (t, x) i— > /ias a continuous modification. 

Additionally, the ODE (|5.4p /ias a unique solution Y f ' x '° G C°([i,T] x 
M m ,IR). Moreover, there exists a constant C > swc/t f/iaf 

| F M,0 _ F M',0| < C | x _ f £ [ 0)T | 5 x y G R m (5 5) 

Proof. Existence and uniqueness of a solution to (|5.3p for e £ (0, 1] follow 
from Theorems 13.91 and 13.101 An argumentation similar to that of the proof 
of (|4.5p yields the uniform bounds in e. The continuity result follows from 
Corollary 14.41 

Concerning Y t,x, ° from (|5.4p . we easily see that g is still a uniformly 
bounded function and F(s,y) = f(s,Xl' x,0 ,y,0) for (s,y) G [i, x 1R is an 
J^-adapted continuous function (in s and y) satisfying (simply use (|4.3p and 
(|4.4p ) a linear growth and Lipschitz condition (in the y- variable). After a 
time reversal, standard ODE theory yields (e.g. Picard-Lindelof's theorem) 
the existence and uniqueness of a continuous solution Y l ' x '® to (|5.4p . One 
then obtains that Y t,x '° is uniformly bounded. 

The estimate for the difference (Y t,x, ° — Y t,x >°) follows easily from the 
same estimate for (X*' x '° —X t,x ,0 ), Gronwall's inequality combined with the 
uniform boundedness of Y t,x '° and (14.41). □ 



The next result states the convergence of {Y l ' x ' s , Z t ' x,£ ) to (Y t,x '°,0) as 
e vanishes. Notice that we can always interpret the deterministic equation 
(g3D as a FBSDE with solution (Y'^O). 



Proposition 5.3. Let Assumptions \4-l\ and 4- 2 hold. Then for any p > 2 



there exists a constant C p > such that 



sup {E[ sup |yW_yt^0| P j +E r/ / |^M^|2 dr yi | < Cp£ P_ 
We emphasize that C p is independent oft, x or e. 

Proof. Let p > 2. We start by interpreting the pair (Y t,x, °, 0) as the solution 
of a deterministic BSDE. This interpretation allows us to use Corollary 13.131 
along with the uniform boundedness of Y t,x,£ and the assumptions to obtain 

||yt,x, B _ y^'Xsp < CpWX 1 '^ - X^Ys^v < C v £P > 

where f3 is related to the uniformly bounded (in e) BMO norm of Z t,x,£ * W 
and the result follows from Lemma 15.11 

From (|3.5p and the results already proved we easily get the sought result 
concerning \\Z^ X ^ - Z^ xfi \\ P Hv = \\Z^^\^. □ 
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The corresponding perturbed PDE problem 

By solving (|5.ip and (|5.3p one is able to solve, at least in the viscosity sense, 
the corresponding PDE, namely in view of (|4.8p and Theorem 14. 101 one has 
for (t, x) G [0, T] x W m 

d t u e (t,x) +C e u e (t,x) + f(t,x,u £ (t,x), y/e{V x u £ a){t,x)) = 0, 

u £ (T,x) = g(x), 

with C e := (b • V)(i, x) + | J2Tj=i( a(jT )i,j(^^ x)d XiXj . And as e vanishes one 
is lead to the first order hyperbolic PDE for (t, x) G [0, T] x M. m 

d t u°(t, x) + (b- Vu°)(t, x) + f(t, x, u°(t, x),0) = 0, u°(T, x) = g{x), 

or after a time inversion, t \— > T — t, with (t, x) £ [0, T] x M m 

d t u°(t,x) = [b- Vu°)(t,x) + f(t,x,u°(t,x),0), u°(0,x) = g(x). 

5.2 The large deviations result 

We now explain the large deviations principle (LDP). For a deeper analysis 
on large deviations we point the reader to |FW98] , |DZ09] , |FK06| and ref- 
erences therein. In this section we extend the results of |DR07| and [Rai06j 
to another class of FBSDE. 

We denote C°([a,b],M. m ) (for < a < b < +oo) the set of continuous 
functions on the interval [0, T] with values in W 71 and we consider in this 
space the uniform norm 

P[a,b}{4>) ■= sup \(j)(t)\. (5.6) 
te[a,b] 

Definition 5.4. The family of processes (Xf)tG[o,T] depending on a parame- 
ter e is said to satisfy a large deviations principle (LDP) with rate function 
S(fy) if the following conditions hold for every Borel set A C C([0,T]) 

limsupelogPpr G A] < inf 

e^0 *£C7(A) 

liminf elogPDT E A] > - inf S(tf), 

£-►0 *e/nt(A) 

where CI (A) denotes the closure of the set A and Int(A) denotes the interior 
of the set A. 

Before we state that X t,x,£ satisfies an LDP, we recall the Cameron- 
Martin space H of absolutely continuous functions. 

Definition 5.5 (Cameron-Martin space). Let fcsN, the space H([t, T], M fc ) 
is defined as the space of continuous functions (j) G C°([t, T], IR fc ) for which 
there exists a square integrable function (ft such that 4> s = <ftt + Jf 0- r dr for 
< t < s < T (i.e. (p is absolutely continuous). 
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We now make a new assumption, stronger than Assumption 14. 1| for the 
result stating that X t,x,e satisfies an LDP. 



Assumption 5.6. Let Assumption \4-l\ hold, b and a satisfy b(t,x) = b(x) 
and a(t, x) = o~(x) for all (t, x) G [0, T] x M m and are uniformly bounded. 

Theorem 5.7 (Theorem 5.6.7 of [DZ09J). Let Assumption \5. 6\ hold. Then 
the solution X t,x,£ of (|5.ip satisfies, as e goes to zero, a large deviations 
principle in the space C°([i, T], M m ) associated to the rate function L x defined 
for any <p G C°([t,T],M. m ) by 

ri r T 

I x (4>) := inf |-y \v r \ 2 dr, v £ H([t,T],R d ) such that 

4> s = x + J b((f) r )dr + J a((p r )v r dr, sG[t,T]|, 
with the convention that inf = +oo. 

Remark 5.8. We point out that if aa T is invertible, then the rate function 
functional can be written in a simplified manner, more precisely 

Ix{4>) . = f \ £ Ql.[<Pr - b(<p r )]dr , if</, € H([t,T],R m ) and <f> t = x 
I +oo , otherwise, 

where Q* is defined by Q^[v] := (v, (aa T )~ 1 (u)v) for all u, v G M m . 

We now define a set of operators from and onto the space of continuous 
functions that, combined with the contraction principle 13 ! (see Theorem 2.4 
in |Var84j or Theorems 4.2.1 and 4.2.21 of |DZ09| ). will allows us to conclude 
that Y l ' x ' e satisfies an LDP. 

Definition 5.9. The operator F e is defined as 

F £ : C°([t,T],R m ) -»• C°([t,T],M) 

0^F £ (0) =n £ (-,0.), 

where for e G (0, 1] the function u e is defined by u 6 (t,x) := Y^ ,x,e the first 
component of the solution to f|5 . 3 j) . And for e = 0, u°(t,x) is given by Y^' x '° 
solution to (15.411 . 



We observe that for all < t < s < T, x G M m and e G [0, 1] we have 
Y t,x,e = F s( X t,x,e^ s ^ = n £( S) xi' x ' ,£ ) (see Theorem IPD. We now state and 
prove the result concerning the LDP satisfied by the law induced by Y t,x ' £ . 



13 In rough, let G and 3 be two complete separable metric spaces and f e : O i— > 3 a 
family of continuous functions in those metric spaces. If lim e _>.o /e = / exists uniformly 
over compact sets, then the contraction principle states that if the valued process X t , 
satisfies an LDP with rate function / then F e = f s (X e ) also satisfies an LDP with the 
same rate and with a raw rate function J(y) = /(/ l_1 '(y)). Moreover, if 7 is a good rate 
function then so is J. See Section 2 of [Var84] or Section 4.2 of |DZ09| . 
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Theorem 5.10. Let Assumption \5.6\ and \4-S\ hold. Then the process Y t,x ' £ 
satisfies an LDP in C°([i, T],M) with rate function I x defined for any ip G 
C°([t,T],M) by 

4(V0 := inf {l x (<fi : <f> G H([t,T],R m ) such that ip = F°(cp)\, 

with the convention that inf = +00 and where I x {-) is the rate function 
from Theorem 15. 71 

Remark 5.11. // the mapping oa T is invertible then the rate function I x 
can be expressed in an alternative fashion using the operator Q* introduced 
in Remark \5.8[ namely for ip G C°([i, T],R) and x G W 71 

T x (ip) = inf j ^ J Q*0 r [<p r - b((f) r )]dr, over 

0GiT([i,T],M m ), <t>t = x, *p s = u\s,<p s ), sG [t,T]}. 

Proof of Theorem 15.101 This proof uses the well known contraction princi- 
ple, see Theorem 2.4 in |Var84| . To apply the contraction principle we need 
only to show that F £ for e G [0, 1] are uniformly continuous operators from 
C°([t, T], W 71 ) onto C°([t, T], M), and moreover, that F 6 converges uniformly 
to F° over all compact sets of C([t,T],R m ) as e vanishes. Notice from Corol- 
lary [4]4]that the mappings x t-t Yg' x ' £ and s 1— > Ys ,x ' £ are continuous (P-a.s.). 

Step 1 - Continuity of F £ : Let e G [0, 1]. From (14~TU|) and f|5~5|) we have 
that x 1 — y u £ (t,x) is Lipschitz continuous. Let KP be a compact subset 
of C°([t, T], W 71 ) and let ((/> n ) ng pj be a sequence of functions from /C° that 
converge to some <p G JC° in the uniform norm, i.e. limn^oo p\t t T](<P n ~ 4>) = 0> 
where p is defined in (|5.6p . From the definition of F £ , (|4.10p . (|5.5p and the 
continuity of K*' we have 

[pi t , T] (F%ct> n )-F £ (<t>))] 2 = SU P |« 6 M?)-« E (*,&)| 2 

s6[t,T] 

= sup |y/'*"' e - r/.^. £ | 2 

s6[t,T] 

<C sup |0«-0 s |2 = C 7^ T] (0"-0)] 2n ^O. 
se[t,r] 

Note that the constant C present in the estimate of Corollary 14.41 is inde- 
pendent of x or (j). 

Step 2 - Convergence over compacts: To prove that F £ converges to F° 
uniformly over the compact sets fc° of C°([t, T], M m ) we start by defining 
the set X := {<p s : (p G JC°, s G [i,T]}, as the set of points consisting of the 
images of the <p G IC in M m for all s G [t,T]. Remark that I is a compact 
set of R m . 
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Now, by Proposition 15.31 (in the particular case of s = t) and taking the 
set I into account, there exists a constant C > such that 



sup 



P[t,T] I 



sup sup |n e (s,0 s ) -u°(s,</> s )| 2 = sup sup |Y/>^> £ - y/.*«.°| 2 
^G/c° se[t,T] <t>eK.° se[t,T] 



<sup sup |y/^ £ -y/' x '°| 2 < Ce 2 ^0 as e ->• 0. 
xgx se[t,T] 

The result now follows from the already mentioned contraction principle. □ 

5.3 PDE perturbation and Burgers' equation with damping 

Here we want to discuss the following PDE for A>0, ofR, i 6R where 

b 



ii:[0,T]xR->R and g G C b ° 



d t u £ + au £ V x u £ + \u £ = eAu £ , u £ (0,x) = g(x). (5.7) 

Or, after a canonical time inversion 1 1->- T — t for t G [0, T] 

<9 4 u £ + eAu e = au £ V x .u £ + Xu £ , u e (T, s) = p(x). (5.8) 

Via the already discussed Feynman-Kac formula we can write (at least for- 
mally) the corresponding FBSDE system for i£R and < t < s < T 

Xl' x ' £ = x + V2e (W a - W t ), (5.9) 

y*.^ £ = ff (X^f' £ ) - / Z^' e dW r - / -== Y^ s Zl' x ' e + Ay*^' 2 dr. 

The appearance of the l/\/2e weight is due to the representation Z = 
(V x ua)(-,X) with cr = V^e; see (02]). This type of FBSDE still fits in 
the setting we discussed in the first sections, but the presence of the term 
1/e creates complications if we consider the FBSDE not just for one single 
fixed e. 

For each fixed e € (0, 1] we obtain from Theorem 14.31 existence and 
uniqueness of the triple (X t,x,£ , Y t,x,£ , Z t,x ' £ ) solving the FBSDE above in 
S p x 5°° x Hbmo for an y P > 2. The real interesting result is the uniform 
boundedness in e of ||y*' x,£ ||s°° . 

Lemma 5.12. Let (t,x) G [0, T] x R then in the above defined framework 

sup ||^* ,:r ' e ||s°° < oo. 

xeK,ee(o,i],*e[o,T] 
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Proof. From the original existence Theorem 13.91 and the lemmata used in its 
proof, we easily obtain for each fixed e £ (0, 1] that 

sup Hl^'^Hs^ < oo. 

x£R,t£[0,T] 

We sketch the argument concerning the uniform boundedness on e since it is 
based on an argument we have already seen in the proof of Lemma 13.161 We 
again use the linearization and measure change argument. Define processes 
et = exp{— At} and bf := aYf /y2e. Since Y £ is a bounded process it is 
clear that b e * W S BMO and we can define a new probability measure Q e 
with Radon-Nikodym density dQ e /dP = £{b £ * W), with relation to which 
W £ = W. — J Q b £ r dr is a Brownian motion. 

Applying Ito's formula to etYf and changing to the measure Q e we obtain 

e s Y^ = e T g(X^ e ) - f e r Z^ £ dW £ 

J s 

Y*> x > s = e- x( - T -^E QS [g(X^ x ' £ )\T s ] 

sup sup || sup |y s ,a: ' e | Loo < sup> I (a?) | < oo, 
xm £e(o,i] se[t,T] xm 

where the last line follows from the uniform boundedness of g and from 
e -A(T-s) < 1 for all s £ [ Q) T\. □ 

This result hints that the boundedness of the solution carries over to the 
limit (as e — > 0). On the other hand, obtaining a similar estimate for Z t,x ' £ 
in the WBMO- norm is impossible without extra assumptions. 

First remark that the limiting PDE of (|5.7p as e — > is the first order 
equation: 

d t v + avV x v + Xv = 0, v(0,x)=g(x), (t, x) e [0, T] x R. 

For this type of transport PDE it is well known that even with smooth 
initial boundary data the solution can develop shocks in finite time which 
translates into an explosion of the spatial derivative in finite time, see |TerlO| . 
An estimate on W x u e (or Z) independent of e would mean that the derivative 
remains stable (non-explosive) as we decrease the regularizing impact of the 
perturbation. 

Getting help from fully coupled FBSDE 

If one assumes that A > K 2 a 2 T then it is possible to prove via fully coupled 
FBSDE arguments that sup £g ( 01 ] ||V X w £ ||oo ^ \/^ K . This condition is, in 
some cases, close to that of Theorem 1 in [TerTO] and it roughly means that 
if a strong enough damping is present in the PDE (via the term \u £ ) then 
the cross term u £ V x u £ does not lead to shocks up to time T. We remark 
though, that |TerlO| works with bounded domains. 
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Even with the mentioned condition these authors have not been able 
to show the result using the decoupled FBSDE theory alone. Notice that 
one could effectively use the theory developed in Proposition 14,61 (with the 
appropriate assumptions) to show that the solution (Y t,x,£ , Z t,x,£ ) (of (|5.9p ) 
is differentiable in x, but the mechanics employed does not lead to estimates 
on the moments of Z t,x,£ or \/ x Y t,x,£ that are independent of e\ 

An alternative FBSDE representation of PDE (|5.8[) follows by interpret- 
ing the term au £ V x u £ not as a component of the BSDE's driver but as a part 
of the second order differential operator, i.e., as the drift of the forward diffu- 
sion (compare FBSDE P~1"]) - (|P]) and PDE (1478]) along with the operator £). 
This leads then to a fully coupled FBSDE, namely take (t, x) G [0, T] x W n 
and s £ [t, T] 



x t,x, £ = x _ a Y*' x ' £ dr + y/2e / ldW, 



ri 



t Jt 

T r T 



Y l s > x < £ = g{X^ e ) - j \Y?> x > e dr - f Z l s ^ £ dW r . 

J s J s 

We give a quick sketch of the arguments for the fully coupled case since 
the theory differs a bit from what has been presented before. Although the 
argumentation is not complicated, a full justification would be lengthy (for 
this already long manuscript); for more on the theory of fully coupled FBSDE 
we point the reader to |MY99| . 

Let us assume that the above coupled FBSDE has a unique solution 
prt^yt^^W) e s p xS°° x Hbmo, and rather trivially (to see) that 
1 1 Y t,x ' £ 1 1 500 and \\Z t ' x ' £ \\y_ BMO are uniformly bounded in their parameters. 
For x, x' £l m ,t£ [0, T] and 5Y = Y f ' x ' £ - Y^'' 5 and 5Z = Z l > x > £ - Z^'' 8 , 
applying Ito's formula to |<5Y" S | 2 for s £ [t,T] and noticing that 

\g(X^ £ ) - giX'/'^l 2 < 2K 2 (\x - x'\ 2 + a 2 (T - t) J* \5Y s \ 2 ds) , 



leads to 



\5Z s \ 2 ds + \5Y t \ 2 



<2K 2 \x-x'\ 2 + 2{K 2 a 2 {T -t)-\) [ \5Y s \ 2 ds-[ 25Y s 5Z s dW s . 

Jt Jt 

Assuming that A > K 2 a 2 T and taking ^-conditional expectation yields 

\8Y t \ 2 = \u £ (t,x) -u £ (t,x')\ 2 < 2K 2 \x-x'\ 2 . 

Surely, this result does not imply differentiability in x, but with some extra 
arguments (mollification of the terminal condition then passing to the limit) 
one easily obtains the spatial differentiability of Y t,x,£ (and hence of u(t, x)). 
From the above estimate, one would obtain |V x it £ | < y/2K uniformly in e. 
This again hints that the estimate would hold in the limit. 
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5.4 Incompressible Navier-Stokes 

While in the previous subsection we were forced in the end to work with fully 
coupled FBSDE in order to obtain the extra result we where searching for, 
we now give an example which simplifies greatly if one uses the decoupled 
FBSDE instead of the fully coupled ones. In this section we look at the works 
by |CS09| ICSlOj where they present a connection between FBSDE and the 
incompressible Navier-Stokes equation. 

More than being a standard PDE, the Navier-Stokes equation describes 
a conservation law (see the equations below) and it is this conservation law 
that is the hard part in the FBSDE representation. We saw in Section U] 
that the FBSDE yields solutions to PDE but the conservation of volume 
condition (i.e. div(u) = {l,V x u) = 0, see (|5.10p ) has to be verified. One 
does not obtain it as a side product of the existence of solution. 

Following [CS09, CS10J we look at the incompressible Navier-Stokes equa- 
tion with a:[0,T]xlMR satisfying for (t, x) G [0, T) x R n 

dtu = — u{l, V x u) + vAii — V x p, div u = 0, (5.10) 
u(0, x) = h(x), div h = 0, 

where v > 0, div(u) := (l,X7 x u) and h G C^(]R n ). The function p denotes 
the pressure field which is assumed to satisfy V x p = K for some constant 
K. Throughout we work with n = 2 but the results presented here are easily 
extended to the general case. 

In view of the calculations carried out in Subsection 15.31 the FBSDE 
corresponding to the PDE (|5.10p (after a time reversal) is for (t, x) £ [0, T] x 

TU2 



Xl> x = x + V2u(W s - W t ), se[t,T] 

Y\* = h(X^) - 2%'dr ~ [ [ K + ^ y ^U, Z 



dr, 



where W = (W l ,W 2 ), X = (X\X 2 ) and Z = (Z l ,Z 2 ). 

From Theorems 13.91 [3~TU1 and |4 . 1 1 easily follows existence and uniqueness 
of a solution (X, Y, Z) to the FBSDE in S p x S°° xHbmo (for any p > 1) as 
well as its connection to the PDE. Moreover since h G Cl, Proposition 14.6 
also holds, yielding the existence of {V X Y^ X , V X Z^ X ) G S p x W for p > 2. 

One of the difficulties in [ CS09| was showing that Y satisfied the free 
divergence condition, namely that for u(t, x) := Y^' x the function u satisfies 
div(ii) = for all (t,x) G [0, T] x M 2 . For the case we present here, this is 
considerably easier than |CS09| . We have 

y x Y^ x = (y x h)(X^)V x X^ x - j\z l r ' x AW r 



f T -L \V X Y^ X (1, Z\: x ) + Y*>*V X ((1, Zp x )) 
Js V 2^ L 



dr, 
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where V X X S = (1, 1). Then, div(ii) = is found by looking at the processes 
U := (l.Va.y* 1 *) and V* := (1, V^Z*-*'*) for i G {1,2}. We have then 



J7 S = (div(/i)) (X 



r 



VrdWr 



T r ii yt,X\ v t,x 



'2v 



2u 



dr. 



Now since div(/i) = and using the linearization and measure change argu- 
ment, i.e. for s G [t,T] define b s := Y s ' x /y/2u, a probability measure Q with 
Radon-Nikodym density £{b* W) and e s := exp{ J/(l, Z*' :r )/v / 2i'dr}. Note 
that since y G 5°° and Z G Hbmo the processes 6 and e are well defined, 
see for instance the arguments used in the proof of Theorem 13.101 
We obtain then 



(IrVrdW® 



U s = Q-a.s. for all s £ [t,T\. 



In particular, Ut = div(Y^' x ) = as we sought. This proves that the solution 
of the FBSDE provides a mild solution to the Navier-Stokes equation (|5.10p . 
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A Proof of Proposition 14.6 

Proof of Proposition \4-6\ The results concerning the forward diffusion are 
standard, see for instance Section 1.2.4 of [DRllJ. 

This proof has 3 steps. We first establish a candidate for the said deriva- 
tive process, then we prove that the candidate is the correct one and we 
finish by showing continuity of the derivatives in the spatial component. 
Throughout let j G {1, • • ■ , m} and t G [0, T]. 

Step 1 - The candidate for the partial derivatives: A quick look at the 
BSDEs composing the lines of (|4.7p shows that the lines are independent of 
each other. So, in view of the assumptions and the already obtained results 
concerning Y x and Z x it easily follows that for each line of the system (|4.7p 
Assumptions A1-A4 of [BC08J hold. Hence we can apply their results to 
conclude that (j4~T|) has a unique solution (V X .Y X , V X] Z X ) G S p x W for all 
p > 2 and j = { 1 , • • • , m} . These processes are the natural candidates for 
the j-th partial derivatives of Y x and Z x . 

Step 2 - The candidates are the correct ones: We now need only to prove 
that the solution of (|4.7p corresponds indeed to the partial derivative of 
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(Y x , Z x ). To this end we need only to show that for any h > (0, 1] such that 
we have for all p > 2: 



lim 



h 



V Xj Y* 



+ 



^x+he 



h 



-v Xj z* 



HP 



0. 



Again we do not give the details, the proof of this result follows from a 
straightforward combination of the mechanics employed in the proof of Corol- 
lary [3J2] and the techniques used to prove Theorem 2.1 of |AIdR07] (alter- 
natively see Section 3 of |DR11] ). 

Step 3 - Total differentiability and the continuity estimate: After having 
established that the partial derivatives exist and satisfy good integrability 
properties, it remains only to argue in favor of total differentiability of the 
map x i — y Y x . This will be done by showing that the partial derivatives are 
indeed continuous at any point x (in the topology of S p x T~L P ). We do this 
by proving the limit result on the difference of gradients. 

This estimate is easy to obtain as it follows from the techniques used 
throughout Section 3. We do not give the full details. Let t G [0, T], x,x' £ 
R m then the BSDE for the difference 5VY = V x Y l ' x - V x Y f ' x ' (we define 
as well 5VZ = V X Z^ X - V x Z t ' x ') is given by 

5VY s = 5r]- f 5VZ r dW r + f [F a + K^5VY r + K?5VZ r ]dr, (A.l) 

J S J S 



with 5rj = 

v*/(.,ef 



K z 



{V x g){X$)V x X% - {V x g){X% )V X X% , K y = V y /(-,6 a 
i and V X Q X ' = (V X X X ', V X Y X ', V X Z X ') we define the process F := 



v x f(-,e x )(v x x* 



v x x x 



+ (v x f(.,e x )-v x f(.,e x ),v x G* 



The results we have seen so far easily imply (see Theorem 14. 3p that K y 
and K z belong to Hbmo with norm uniformly bounded in x. From this 
structure it is easy to see that the technique used to prove Corollary 13.131 
can be applied to this situation, yielding for a (3 > 1 related to the BMO 
norms of the relevant processes and a constant C independent of x,x' that 



\SVY\\ P SP < CE 



N p/3 + 



\F s \ds 



P j3- 



We point out that although (lA.ip is a multidimensional equation and our 
trick consists of a measure change, the measure change is uniform across all 
entries of VY", namely for each entry it is always the process K z * W that 
is at the origin of the measure change. 

Given the already proved (see Corollary I4.4[) continuity of X x , Y x and Z x 
in x, the continuity of the derivatives of / (by assumption), one can easily 
conclude via the dominated convergence theorem the sought result. This 
type of arguments have been used in Section 3.1 of |DR11| or [IDR10J. □ 
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B Proof of Theorem 14.101 

Proof of Theorem\4-l(J\ Let (t,x) G [0, T] x R m . We have seen in previous 



results that the process Y t,x has continuous path as well as a continuous 
dependency in x (see Corollary I4.4[) . Uniform boundedness of u follows from 
that of Y implied by Proposition 14,31 

To prove that u is a viscosity solution we need to prove that it is a 
viscosity sub- and supersolution. We prove only that it is a subsolution 
since the supersolution part of the proof if very similar. 

To prove that u is a viscosity subsolution, let tp G C 1 ' 2 ([0,T] x M m ,R m ) 
and take (t, x) G [0, T] x M m such that (t, x) is a local maximum of (u — tp). 
We assume that without loss of generality that they touch, i.e. that u(t, x) = 
ip(i,x). Suppose thaif^l 

d t <p(i,x) + Cip(t,x) + f(t,x,u(t,x), (V(pa)(i,x)) < 0, 

and we now argue by contradiction. 

Since [0, T] x R m is a separable space we can define a neighborhood of 
(t, x) in the following way: let < a < T — t be such that for allt < s < i+a, 
\y-x\< a, u(s, y) < cp(s, y) and 

d t ip(s, y) + &p{s, y) + f(s, y, u{s, y), (Vipa)(s, y)) < 0, (B.l) 

and define 

r = inf{s > i : \X 1 / - x\ > a} A (i + a). 

The pair 

(Y s , Z s ) := (Y*&, l [t>] (s)zl^) , s G [t, t + a] 
clearly solves the BSDE on a] 

Y s = u(t,X^)+ [ t+0 l [i}T] (r)f(r,X^,u(r,X^),Z r )dr- Z r dW r . 

J s J s 

On the other hand, since if G C 1 ' 2 we have by Ito's formula that 

(Y S ,Z S ) := (<p(8,xl* T , l %T] (s)(Vipa)(s,X%)), s G a] 
solves the BSDE for s G [t, i + a] 

Y s = <p(r, X l f) - [ t+0 l[t, T ] (r) {d t y + C<p) (r, X^)dr - Z r dW r . 

J s J s 



1 Compare with the definition of viscosity subsolution: Take u £ C([0, T] x I 
tp G C 1,2 ([0,T] x R d ) and (t,x) e [0, T] x R m . Then it is a viscosity subsolution of (g] 
if for all x it holds that u(T, x) < g(x) and if for an <p and a pair (t,x) which is a local 
maximum of u — <j>, it holds that —dtf(t, x) — Cipit, x) — f(t, x, u(t, x), (V<po~)(t, x)j < 0. 
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At this point we need only to use a comparison result that allows us to con- 
clude that Y < Y and hence that u(t, x) < (p(t, x) to find the contradiction 
with our assumption that u(t, x) = (p(t,x). Unfortunately we cannot apply 
directly our comparison Theorem I3.12[ because we do not know if the BSDE 
for (Y, Z) satisfies its assumptions. 

We start by remarking that since ip E C 1 ' 2 and within the neighborhood 
around (t, x) defined above, tp(r, Xr' x ) and (dt(p+C<p)(r, Xr' x ) are continuous 
functions on a compact set and hence bounded! This allows us to trivially 
conclude (e.g. Theorem 13. 7p that (Y, Z) E 5°° x Hbmo on the time interval 
[i, i+ a]. It is also clear that (Y, Z) belongs to 5°° x %bmo simply because 
it is the solution of ([4TT]> . (jl^l . The BSDE for the difference Y - Y is 

Y s -Y s = u(t, X l f) - tp( T , X^) - Z r - Z r dW r 

J s 

+ f t+a i [lr] (r)[f(r,X^,u(r,X^),Z r ) + (dt<p + £<p) (r, X*>)] dr. 

J s 

Adding and subtracting a term f(r, Xl' x , u(r, Xr ,x ), Z r ) to the driver of the 
BSDE for the difference Y — Y we can apply arguments similar to those 
of the proof of Theorem 13.121 and conclude from the fact that u < ip and 
the strictness of (IB. 1|) that Y < Y, i.e. u(t,x) < (p(t,x). This contradicts 
the assumption and hence we can conclude that u(t, x) = Y t t,x is indeed a 
viscosity subsolution of the PDE. 

The arguments to show that u is a supersolution are similar. A combi- 
nation of the two resutls yields that u is a viscosity solution of (j4.8p . 1 □ 
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